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ABSTRACT 


This thesis presents an interactive naval radar sodel 


which computes radar detection in the presence of land 


masses, using a parametric terrain descripticn. 
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I. ZNTROPUCTION 


"This thesis presents a novel naval radar nodel which 


computes radar detection in the presence of land masses. The 


model is an interactive computer program which accepts 


scenarios and radar parameters from the user and displays a 


map of the area indicating where targets can and cannot be 


detected. The resulting map can be displayed at the user's 


computer terminal or printed offline. 


Major capabilities of the model are: 


parametric terrain description 


user friendly input and output 


beach return masking 


radar shadowing 


side lobe aasking 


The program is written in FORTRAN IV H (extended), to be 


executed on an IBM 3033 using an IBM 3278.2 video computer 


terminal. A data file which contains a parametric terrain 


description aust be prepared before using this program, but 


all other required input is prompted at the terminal. 
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II. MOTIVATION 

The Warfare Environment Simulator (WES) used by Command, 
Communications, and Ccntrol personnel is a large scale con- 
puterized naval wargame. The land dispiayed by the systen 
does not affect the radar detection probabilities. This is 
one of the artificialities in the game. 

During Exercise BRIGHT HORIZON '81 the tactic of con- 
cealing small vessels in fjords and among islands to prevent 
radar detection until a surprise attack could be made was 
demonstrated to be effective. This technique, while it 
could be anticipated, can not currently be modeled in WES. 

North Atlantic Treaty Organization naval warfare scenar- 
ios virtually always include the proximity of land. United 
States Navy scenarios are primarily open ocean. It seems 
necessary to incorpora*+é NATO probleas into USN procedures. 
With this in mind, a method of encoding terrain data was 
combined with a radar model to produce a method of producing 
maps of the area in question with the concealed areas 


displayed. 
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III. MODEL 


A. WARFARE ENVIRONMENT SIMULATOR 
1. Basic system 
WES is a computer program that siaulates a naval 
warfare environment in a large computer system, enabling 
personnel to engage in realistic wartime scenarios without 
the expense of actually losing ships and lives {[Ref. 1]. 
2. Radar Equation 
The WES model uses the standard radar equation nodi- 
fied to handle ducting and clutter effects. Equation (1) is 


used to calculate signal excess. 


SE=P* #¢2G42W+TCS-4DR-B~NF-L-C (1) 
SE = signal excess, dB 


Pt = peak transmitted power, dB//w 
G = antenna gain, dB 

W = wavelength, dB//ca 

TCS= target cross section, dB//nm 

D = ducting factcr 

R = target range, dB//na 

B = receiver noise figure, dB 


NF = noise factor, 45 
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clutter factor, dB 
L = system loss factor, dB 

System loss factor includes antenna pattern loss and 
atmospheric absorption loss. Antenna pattern loss takes 
into account the change in illuminating energy lavels 
brought abou by the lobe shape of the radar main beam pat- 
tern. The internal system losses are assigned to he 1.5 dB, 
a typical value for macst radars. Atmospheric absorption loss 
is added to this. Atmospheric absorption loss is frequency 
dependent, and is assigned as 1 d@B around 300 MHz and 3.5 dB 
at 5000 MHz, at target ranges typical to naval radars. 

The clutter factor describes the losses due to sea 
clutter. It is taken to be 10log(Hw), where Hw is the signi- 
ficant wave height in feet. 

The ducting factor is used to describe the effects 
of ducted propagation, especiaily surface evaporation ducts. 
An approximate fit of the IREPS model of a leaky wave guide 
is used.! A strength of zero indicates no ducting, anda 
radar horizon limitation is evoked. If a duct is present, 


this restriction is lifted. A strength of five corresponds 


&Integrated Refractive Effects Prediction Systea 
developed by Naval Ocean Systems Center and ap emented on 
the Hewlett-Packard model 9845 desktop. calculator. It is a 
shipboard anys cgaentad data processor and dis lay systen 
designed to aid in the assessment of the iapa ofl 
atmoSpheric refraction effects on Naval EM systens. 


wer 
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to a "perfect", lossless duct with free space propagation 
within its boundaries. Values for D are selected from Table 


(1). 


Table 1; Ducting Factor 


rgngt ht jp ay 8 dg ig 


st 

The noise factor is used to include jaasing. It is 
not used in this sodel. 

See subsection 3 below for the aethod of processing 
signal excess. The method used in determining the radar tar- 
get cross section of ships is discussed in Chapter III, sec- 
tion c, subsection 3. 

Por this model an additional term is included. The 
“obstacle gain" is calculated and subtracted from the right 
hand side of the equation.2 

3. Signal Excess aode} 

Signal excess is converted to probability of detec- 
<ton on the basis of false alarm number and pulse integra- 
tion as discussed below. 

False alarm number is the numoer of radar signal 
pulses that are received before a non-sigaal aoise pulse is 


4 weer obstacle gain see chapter III, section c, subsec- 
on 6. 
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received. WES, as well as this aodel, offers the use of four 
different values. 

The number of pulses illuminating the target is 
evaluated by use of eguation (2). 


60 (PRR) (BW) 


NP Sa) (ARR) (2) 
where: 
Np = number of pulses 
PRR = pulse repetition rate, pulses per second 
BW = horizontal beam width, in degrees 
SW = angular width of swepth sector, in degrees 


ARR = scan rate, in scans per minute 
The distribution of signal excess can be closely 
approximated by a normal distribution with a mean defined by 
the number of pulses integrated, probability of false alara, 
and a standard deviation of 7458. A large takle contains the 
mean of the signal excess distribution (mu) in terms of nus- 
ber of pulses integrated and probability of false alarn. 
Conversion of computed signal excess (SE) to prob- 
ability of detection is accomplished in WES as follows. 
The signal excess is normalized by equation (3). 
z 2 (se-mu)/sd (3) 


The z is used to enter a table of cumulative standardized 
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normal values to gat the probability of detection which 
is used in the WES siaulation. 

This model runs the process in reverse: given the 
desired probability of detection, the number of pulses, and 
an assuaed standard deviation the standardized normal prob- 
ability distribution is used to find a corresponding signal 
excess, which can then be compared to the realized signal 
excess to determine if the target will be detected within 
the constraints given. 

4. Conditions and Assusptions 

A number of constants are established in the WES 
model, as well as in this modified version. Many of then 
will be discussed in the sections in which they are singu- 
larly appropriate, but some of the more general ones will be 
mentioned here. 

A 4/3 earth radius is used in calculating the dis- 
tance to the radar horizon to coapensate for standard ataos- 
pheric refraction. 

Receiver noise figure is taken as a constant value 
of 5.5 dB, a value which is prigarily thermal noise under 
standard conditions. 

The return froma land is not directly modeled. The 


geometrical cross section of a portion of parametrized 
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terrain is assumed to be the radar cross section. Data is 
unavailable on the reflectivity of varying types of surface 
to compute the true radar cross section given the geoaetric 
cross section, so a reflectivity of 1 is assumed. A slope 
of at least .02 is required for beach return, when checking 
from the seaward side, and a slope of at least .3 is 
required for return in the sidelobes. In either case, the 
land must be at least as high as the target vessel. These 
values were determined by experimenting with the geometrical 
cross section of a terrain sample. The steeper slope 
required for the sidelobe masking is caused by the lower 


gain present in the sidelobe. 


B. SNOLER-MELLS MODEL 
1. Basic System 

In 1979 Josef Saoler wrote his Master's thesis on an 
"Operational Lanchester-Type Model of Small Unit Land Coa- 
bat" (Ref. 2]. It is a time sequenced, deterainistic, batta- 
lion-level, force-on-force model implemented on a digital 
computer. It contained a method cf modeling terrain devel- 
oped by Christopher James Needles in March, 1976 for his 
Master's Thesis "Parameterization of Terrain in Army Combat 
Analysis" [Ref. 3]. In September 1980 Glenn &. Hills 


atteapted to overcome shortfalls in the original model, and 
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added several enrichments to provide added user flexibility 
(Ref. 4j. A user*s manual is provided for this model on a 
permanent disk in the W. R. Church Computer Center. These 
Fortran programs conduct their Army battles over land, using 
Needles‘ terrain model and Professor Hartman's elevation and 
line-of-sight routines.3 These routines were used to saodel 
the land in this naval model. Only minor modifications were 
required, partially tc adjust to the change in scale and 
partially to remove some strictly land effects. 
2. Terpain Model 

In 1976 Christopher James Needles presented and 
evaluated a methodology of parameterizing terrain for use in 
land combat analysis (Ref. 3J. This was a shift from the 
traditional method of digitizing data compiled from terrain 
and interpolating. He described a method by which terrain 
could be created mathematically by using a aodified bivari- 
ate normal probability density function. 

The common forma of the bivariate normal density 


function is given by equation (4). 


3see sections three and four of this chapter for ELBEV 
and LOS. 


eae. 
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The equation has been modified by making the noraal- 
izing constant equal the the naxinua elevation of the 
desired terrain. The resulting equation has sufficient par- 
ameters to model many different types of “hills", and by 
carefully combining these hills different types of terrain 
can be modeled. 

This method of terrain modeling is currently being 
used by the STAR model written in SIMSCRIPT II.5, but the 
Fortran version of the method has also been used, in the 
Smoler-Mills model previously described. *¢ 

See Figure (1) for an illustration on how the hills 
are constructed. 

3. Elevation Routine 

The purpose of the function ELEY is to coapute ter- 
rain elevation for given X, Y coordinates. 

The elvation routine was initially provided to Josef 
Saoler by Professor James K. Hartman for use in his land 
combat mod3l, and subsequently by Millis in his modification 
cf Smoler's model (Ref. 5). 

Tha point defined py the coordinates is exanzined to 


determine if any of the "hills" are present, and if so which 


“Simulation of Tactical Alternative Response, a 
aaa combat agaodel developed at the Naval Postgraduate 
chool. 
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Figure 1: Terrain Structure 


ones. Only those which are effective are avaluated, siniaiz- 


ing tha computer time required for each point. Given the 


2c 
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index of the “effective hill", the appropriate paramenters 
are selected and applied to the previously defined rnodified 
bivariate normal distribution to derive the terrain eleva- 
tion at that point. 

In a second elevation routine, ELEVG, the gradient 


components of the equation are computed as well as the ele- 


ne 2 ey ear ee 


vation, allowing the slope to be determined easily. 

The procedure to add tree height to the top of the 
hills to increase elevation has been reaoved for this radar 
model. 

4. Line of Sight Routine 
The line-of-sight routine was also developed by 

Professor Hartman (Ref. 5]. Its development and use follow 

| 


Closely those of the elevation routines. 


Line of sight (LOS) is a purely geometric computa- 
tion assuming perfect visibility. The result of the LOS coa- 
putation was initially the percent of the vertical height of 
the target visable to the searcher, and the percent visible 
of the searcher by the target. In addition, the range and 
elevacion of the point which causes the target to be totaily 
ebstracted is now returned. 


The basic procedure is: "Pind the lowest sight line 


Extend this line to the 


froa the searcher over the terrain. 


ee eee, 


target's position and compare its extrapolated height to the 
target's elevation. Thus compute the percent visible." 

In as much as the searcher and the target are both 
ships, the searcher is assumed to be at sea level. If duct- 
ing is present, the target is also assumed to be at sea 
level. If ducting is not present, an additional routine coa- 
putes the amounts by which both the obstructing terrain and 
the target are lowered by the earth's curvature. 


As in ELEV, forest calculations have been removed. 


C. AUGMENTATION 

1. Introduction 

i‘ The STAR LOS and ELEV routines have been combined 
with the WES radar equation to get a routine to calculate 
detection behind and next to land. 

A separate routine was written to calculate beach 
effects, side lobe effects, and a special routine to fro- 
vide output on an IBM 3278 computer video terainal. 

After accepting preliminary terrain, radar, target, 
and anvironmental data the prograa samples sites a desig-~ 
nated distance apart from a specified reference point to 
create a matrix containing either elevation or detection 


data, to be examined at the terminal, with an option to 


print it on a line printer. 
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Pigure 2: Flow Chart 


Figure (2) presents a flow chart of the overall 

program. 
2. Input Routines 

The following subroutines were developed to input 
the information required by the programs. 

a. Terrain Data 

A number cf different items ate required to 

define the terrain. The elevaticn of the reference plane 
must be entered, along with the total number of "hills. For 


each hill, the X and Y¥ coordinates, its saxinum elevation, 


the standard deviation in both the X and Y direction, and 
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the correlation between them must be entered for each hill. 
This defines the hill, in as much as these are the appropri- 
ate parameters in the bivariate normal distribution 
described in the secticn on the terrain nmodel.$ 

The data for the terrain is assumed to be con- 
tained in a data file accessable by the virtual machine. Two 
routines adapted from ones writted by Mark L. Yount perforsa 
the necessary file definitions, after requesting filename, 
filetype, and filemode (filemode defaults to a) at the 
terminal. 

Information for up to 100 hills covering over an 
area 200nm by 200nm may be entered this way. 

b. Radar Data 

Radar information may be entered in either of 
two ways: First, it may be stored in a data file, and loaded 
similiar to the method utilized to load terrain data. Sec- 
ond, it may be entered at the terminal, in response to ques- 
tions on the screen. Table (2) lists the information 


required. 


SPor BYN equation used in the terrain model, see chap- 
ter III, section 8, subsection 2. 
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Table 2: Required Radar Data 


| 
| NAME OF RADAR SYSTEM 
i PEAK TRANSMITTED POWER * aay 
RECEIVER IP.3ANCDWIDT 
PULSE REPITION RATE at 
HORIZONTAL BEAM WIDT E eauce 
VERTICAL, SEAM WIDTH A2z6 Eh 
ANGULAR WIDTH OF SWEPTH eecr R (DEGREES) 
SCAN RATE (oP) 
PULSE LENGTH s ICROSECCNDS) 
FREQUENCY {#8 TZ) 
ANTENNA HEIGHT (FEE } 
S-E COORDINATES OF RADAR (NM) 
i NUMBER OF LOBES (UP TO 
STRONGEST # LOBES (DEGREZS CW FROM MAIN LOBE) 
‘ (notes the *#" above represents the number of lobes desired) 
c. Target Data 
The target data is entered in the same subrout- 
, ine. Table (3) lists the required information. 
| 


Table 3: Target Data 


IDENTIFIZR FOR TARGET 
VERTICAL SIZE OF TARGET feet 
TARGET MAXIMUM DISPLACEMENT (KILOTONS) 
ad. Other 


Additional information is required for the pro- 


gram. It consists of environmental and probabilistic data. 


Tabia 4: Environnental and Probabilistic Data 


DUCTING STRENGTH a ag FOR NO DUCT TO S FOR PERFECT DUCT) 
PROBABILITY OF FALSE ALA 

REQUIRED PROBABILITY OF Den ect ron FOR SINGLE SWEEP 
SIGNIFICANT WAVE HEIGHT (FEET) 


Table (4) lists these additional requirements. 


A copy of initial conditions is then available 
to be printed on the line printer upon request. Figure (3) 


is an example of this listing. 
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Finally, information relevant to the actual 


plotting is requested. Simply enter the nuaber of nautical 

miles between sample points desired and a reference point on 

the overall map and the routine commences its calculations. 
3. Computing Target Cross Section 

The Target Characterization Branch of the Naval 
Research Laboratory has aade carefully controlled seasure- 
ments of the radar cross sections of a number of naval 
ships. A simple empirical expression has been obtained (see 
equation (5)) that expresses the cross section as a function 
of the displacement and the frequency. 

tes=52f*at (5) 
Where: 
tcs = Taryet cross section is square aeters 
£ = radar frequency in megahertz 
a = displacement in kilotons (full load) 

This expressicn provides reasonable approximations 
for the sicrowave band common to most naval radars, and over 
the displacement range of 2000 to 17000 tons {Ref. 6]. 

4. Beach Retusa sJasking 

Each segment of land (elevation greater than 0.5 

meters) is checked to determine if it is a beach. This is 


done in subroutine BEACH by checking every adjacent saaple 
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point in the grid and determining if any of thea is ocear 


(elavation less than 0.5 meters). If this is the case, it is 


assumed that a potential target can get next to it. Next 


Subroutine BRIN estimates the geometric cross section of a 


section of beach. This is done by first determining the area 


within one resolution cell. The size of the cell is deter- 


mined by taking the area between two circles at the given 


range separated by the range resolution and between two 


tadial lines separated by the angular resolution. Subroutine 


ELEVG is then called to determine the slope of the section 


of beach at the point of concern as seen by the radar. The 


effective area of the resolution cell is determined by tak- 


ing the sine of this slope and multiplying it by the resolu- 


tion cell total area. See Figure (4) for an illustration. 


If the effective cross section determined in the 


above manner is greater than that of the target vessel, it 


is assumed that any return from the beach would aask the 


return from the target. 
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5. Main and Side Lobe Interference 


LAND is a subroutine which checks the main lobe and 


given major side lobes for the presence of land. A proce- 


dure similiar to BEACH is used to determine the proxigity of 


rr eee eneny 


A is range resolution 
Bis angular resolutjon 
C is angular resoiutiocn 


| 
| 
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| TOP VIEW SIDE VIEW 
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land, and once *his is established ELEVG is used to check if 


Figure 4: Geometric Area 


there is sufficient slope for ctadar return. 

Sid2 lobes are checked by using ELEVG to check for 
land at the same range as the point to be evaluated in the 
direction of a major lobe when the target is centered in the 
side lobe. 

The presence of land of sufficient height and slope 
in any of the inspected locations indicates that masking 


takes place. 


6. Radja Shadowing 
After calling the line-of-sight routine, the prograa 
*reats the contact in one of three ways: 
a. Case 1 
There is no obstruction. The obstruction gain is 
set +0 zero. 
bk. Case 2 
There is only partial obstruction. The target is 
at the interface between shadow and illumination. The obsta- 
cle gain is set at «he log of an eighth, in accordance with 
+*he basic optical interference technique. [Ref. 7] 
c. Case 3 
The target is totally obstructed, in which case 
subroutine OBGAIN is called to determine obstacle gain. 
OBGAIN calculates the "obstacle gain" created by 


the cadio waves being diffracted over the hill tops. A 


"knife edge" obstacle is assumed, and the calculations are 


performed using Picquenard's radio wave propagation equa-~- 
sions fRef. 3]. If the hill is not the highest of the 
radar, <he hill, and the target it is assumed that the tar- 
get is again at the interface and the same value mentioned 


earlier is returned, otherwise the calculations continue. 


In any case, the obstacle gain is treated as an 


values in the 


attenuation and subtracted from the other 


equation used to calculate signal excess. 


7. Display Routine 


RMAP is a version of TMAP modified to allow the use 


of special symbols. TMAP has been submitted for the NONIMSL 


LIBRARY at NPS. The routine interactively takes an array and 


prepares a contour map for the vicinity of a given point. 


The output is on a terminal or on the printer, at the option 


of the user, interactively taking input from the keyboard. 


Passed arguments: A is the array to be plotted 


N and M ara its dimension 


SFNM is the sumber of nautical niles 


between points. 


Other arguments are requested at the terminal. the 


program is a modified version of Professor Gilles Cantin's 


program CTRMAP version of OCT 1, 1969, adapted for terminal 


use. 


Figure (5) is given on the screen immediately before 


comgencing the plots, either initially at the terminal or 


printed at the top of the paper printout. 
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Syabology for Plot 


Figure 5: 
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Note: On Figure(6), the “$" representing the radar's position is off the screen. 
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Figure (6) is an example of what appears on the ter- 


minal. The paper listing gives the entire aap area, with 


symbology as before. 


PLOT CENTERED ON 45.00NM SOUTH AND 45.00NM EAST 


J 
I B *eekeres 
D Jd UW WU J D + &eeaeeaaenne 
*#C F S W WUS HF CBA*##k424484492 
* E QSU WW USQ E B *#eeeeeeeeugankee 
**® 8FGI TU UT E CRERKESSERKASKSRERLEER 


Pigure 6; Screen Display 


In the screen display given, the points are one nautical 


mile apart, with the radar located approximately 20nm to the 


northwest of the center of the island. The letters roughly 


centered in the screen correspond to terrain elevation, as 


coded in Figure (2). The '** on the left edge correspond to 


a steep beach, where the beach return is stronger than the 


raflection from the target. The '** to the right define the 


area shadowed by.the island. The clear area around it is sea 


level. 


8. Rerun with Charged Data 


This rcutine allows the model to be rerun after 
changing a selact number cf the parameters. Figure (7) is 


the screen display which is the menu of allowable changes. 


Siaply enter the number corresponding to the desired change 


and then the new value. The routine will offer to print the 
new parameters and then construct another matrix containing 


elevation and detection data for exagination. 
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IV. #UTGRE 


A. POSSIBLE UTILIZATION CF MODEL 

This model has applications in planning and training for 
Naval Operations. 

For training purposes, this model could be incorporated 
into WES or other sea combat models to better reflect detec- 
tion in areas where land is present. For such uses, the 
input and display routines could be trimmed, and only the 
point in question would be analyzed, uot an entire map 
array. 

For planning purposes, ship courses could be selected to 
provide maximum radar protection from land shadowing, or the 
inverse, to allow minimum terrain iaterference for stationa- 
ing radar guard ships, could be found after parameterizing 


the applicable areas. 


B. FUTURE ENHANCEMENTS 

Parabolic cylinders could be utilized to aodel tha 
refcaction effects instead of knize-edges to improve the 
model somewhat (Ref. 9]. 

Tha input and output routines could be aodified to take 
better advantage of facilities available (plotting 


terminals, for instance). 


A method to find the beach, instead of identifying if a 


particular point is within a given distance of the sea, 


would enable better ccastal coverage. 

Radar return is often observed behind obstacles under 
conditions in which simple knife-edge refraction cannot be 
responsible. The causes for these effects should be found 
and incorporated into the model. 

For use in other gwodels, different subroutines may be 
adapted for the particular model being constructed. Some 
subroutines will have to pe entirely rewritten, especially 
those involving beach return which use mapped data already 


calculated, to handle single-point problems. 


The following is an entire listing of tne radar nodel. 
It has been specially written to be run on the I5M 3033 on a 


IBM 3278.2 video terminal, the facilities available at the 


Naval Postgraduate School. Appendix B contains and explana- 


tion of each of the subroutines. 
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APPENDIZ B 


SUBROUTINE DIRECTORY 


The following subroutines are used in this project: 


A. BEACH 

This subroutine determines if a point of land is adjacent to 
the sea. 

IX and IY are ¢he X and Y coordinates in array A 

A is the map array, dimensioned N by & 


KOAST is the returned value 


B. BRIN 
This routine calculates the geometric cross section of an 
area of land the size of a resolution cell which may contain 
a target. 

X1,Y¥1 are the position of the radar 

X2,Y2 are the positicn of the beach 

RRES is the range resolution 

ARES is the angular resolution 

TMAC is the terrain elevation 


ECS is the effective radar cross section of the beach 


C. ECURVE 
This subroutine computes the elevation reduction caused by 
the earth's curvature present at a given range. 

BER is the effective earth radius in meters 

RANGE is the distance traveled along the earta 

DROP is the loss in elevation resulting from the earth's 


curvature 


D. ELEV 
This routine determines the elevation at a given location. 
X,Y are the coordinates of the point 


TMAC is the returned elevation 


Ee. ELEVG 
This reutine is similiar to ELEV, but it also calculates the 
gradient components. 


GX AND GY are the calculated components of the gradient. 


F. GETSE 
This routine processes radar data to compute the signal 
excass returned to the radar from a contact under the speci- 


fied conditions. 


SE is the returned signal excess 


G. INPOT 
This routine reads in the radar, target, environmental, and 
detection data, either from the terminal or from a data 
file. (the variables are listed in order of appearance in 
the subroutine) 
RDRS$ name of the radar syste 

peak transmitted power 

receiver bandwidth 

pulse repetition rate 

horizontal beam width 

vertical beam width 

angular resolution 

antenna gain 

angular width of swept sector 

scan rate 

pulse length 

frequency 

radar antenna altitude 

radar anterna altitude 

coordinates of radar 


wavelength 


number of pulses illuainating a target per sweep 


identifier for target 


vertical size of target 


vertical size of target 


target maximua displacement 


target radar cross section 


ducting strength code letter 


probablity of failure code letter 


reguired probability of detection for a single 


sweep 


signal excess required to give a desired PDET 


significant wave height 


clutter factor 


| RL system loss factor 


RNP receiver noise figure 


| RMIN miniaum range 


RRES range resolution 


RMAX maximum unambigious range 


RHORIZ radar horizon 


H. INTRO 


This is the initial page display. 


I. KOVER 


This routine is only used internally by LOS 


J. LAND 
This subroutine determines if a point would be masked by 
return from land in either the main beam adjacent «to the 
given point, or by land at the appropriate point in a side 
lobe. 

MASK is the returned value, one if there is masking, zero 


otherwise. 


K. LOS 
This routine calculates the line-of-sight in terms of a 
fraction visible for observer-target pairs. (the variaples 
are listed in the order of appearance in the subroutine) 
XA,YA (X3,Y3) are the X,Y coordinates or the field for a 

and 3. 

(TMACB) are terrain elevation for A (B) 

(TMICB) elevation of A (B). 

(SIZEB) vertical dimension of A (B). 


(LBTOA) indicator variable for LOS calls. 


HHW is the height of the obstructing hill (if any). 


WRANGE is the range from the radar to the obstructing hill 


(if any). 


VISPRA (VISPRB) fraction of SIZEA(SIZEB) which can be seéen 


by B(A). 


OBGAIN 


L. 


This routine calculates the “obstacle gain created by the 


radio waves being defracted over hill tops. 


elevation of the radar 


H1 


HO elevation of the hill 


H2 elevation of the target 


D1 range from the radar to the hill 


D2 range from the hill to the target 


RLAMDA wavelength 


OBGN returned obstacle gain 


M. $PARS 


It splits a double-precision 


This routine is by Mark Yount. 


string down into its separate words. 


N. RANGE 


This routine finds the range between (X1,¥1) and (X2,¥2) in 


a plane. 


RANGE is raturned. 


O. RMAP 

This routine prepares contour maps with special symbols to 
be displayed either on the terminal, in response to loca- 
tions keyed in, or on a line printer. RLAP is composed of 
three subroutines; PREMAP, PSTMAP, and PRTMAP. 

A is the array to be plotted, dimensioned N by 4M 
SFNM is the number of nautical miles between sample 


points on the grid. 


P. SELECT 
This routine defines the data files so that the program may 


access them. Written with assistance from Mark L. Yount. 


Q. SETUP 

This routine treads in the terrain data from a predesignated 
data file. (the variables are listed in the order they 
appear in the subroutine) 

NAILLS total number of hills 

BASE surface reference plane 

XC(I) and YC(I) are the x,y coordinates of the center of 
the aill 

PEAK(I) is the maxiaum elevation of nill I. 

SX(I) and SY¥(I) are the standard deviations in X and Y of 


hill [. 


719 


| 
| 


RHO (TI) 
LST 


NHL 


NHTOT 


LISTH 


correlation between X and Y for hill I 

pointer on where to start reading on LISTH 

array holding how many hills are in corresponding 
grid square 

total number of hills in LISTH 


a list of which hills affect which gridsquare 
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